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PREFACE 


This  report  describes  the  development  and  application  of 
a method  of  analysis  suitable  for  preliminary  design  of  earth 
penetrators.  Only  normal  impact  events  are  considered.  The  work 
was  performed  by  the  School  of  Aerospace  Engineering,  Georgia 
Institute  of  Technology.  Professor  L.  W.  Renfield  served  as 
Project  Director  and  Co-Principal  Investigator.  Professor  S.  V. 
Hanagud  was  the  other  Co-Principal  Investigator. 

The  Project  Officer  was  Major  David  R.  Spangler,  SPSS,  Defense 
Nuclear  Agency. 
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INTRODUCTION 


Overview 

The  work  described  herein  is  part  of  a larger  program  sponsored 
by  the  Defense  Nuclear  Agency  to  develop  the  base  for  earth  penetration 
technology.  The  objectives  of  the  overall  program  are  to  assess  the 
state-of-the-art  in  design  analysis  prediction  techniques  for  realistic 
earth  penetrating  vehicles  (EPV's)  and  to  establish  a groundwork  of 
experimental  data  obtained  under  controlled  conditions.  Several  parallel 
analytical  prediction  studies  have  been  undertaken,  along  with  two  well- 
instrumented  penetration  tests. 

The  tests  have  been  performed  by  the  Avco  Corporation  in  so-called 

★ 

reverse  ballistic  fashion.  Two  normal  impact  events  were  staged  for 

sandstone  target  material  and  approximately  1800  feet  per  second  impact 

velocity.  The  penetrators  used  in  the  tests  were  designed  and  built 

by  the  Sandia  Corporation^  for  reverse  ballistic  testing;  they  are  half 

** 

scale  and  shown  in  Figure  1.  The  penetrators  were  thoroughly  instrumented 

with  36  strain  gages  and  3 accelerometers  installed  at  locations  that 

2 

were  both  important  and  accessible.  The  data  taken  is  currently  avail- 
able in  preliminary  form  as  Reference  3. 

A simplified  method  of  analysis  for  use  in  preliminary  design  of 
EPV's  is  described  herein.  In  addition,  strain  and  time  history  predic- 
tions for  key  locations  on  the  primary  structure  obtained  by  this  method  are 


★ 

This  type  of  test  is  conducted  with  a stationary  penetrator  and 
a moving  target  block. 

★ ★ 

In  the  sequel,  we  will  refer  to  these  simply  as  the  "penetrators". 
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The  predictions 


are  presented  for  comparison  with  the  penetrator  tests, 
are  obtained  from  simple  closed  form  analytical  expressions  and  are 
in  excellent  agreement  with  available  experimental  data.  This  indicates 
that  the  method  of  analysis  is  both  simple  and  effective. 


Scope 

The  intent  of  our  work  has  been  the  development  of  a simplified 
method  of  analysis  for  use  in  preliminary  design  of  the  penetrator 
primary  structure.  Our  objectives  have  been  simplicity,  rapid  turnaround 
and  low  cost,  in  addition  to  sufficient  accuracy  for  adequate  prediction 
of  primary  physical  processes  and  associated  variables.  Our  major  con- 
cerns are  stresses  and  accelerations  in  the  vehicle  body  itself  and 
stress  transmission  to  the  payload  shell  or  liner.  (Refer  to  Figure  1.) 

Structural  dynamic  analysis  methods  for  EPV's  presently  in  use 
employ  large-scale  numerical  analysis  computer  codes  based  upon  finite 
difference  or  finite  element  structural  models  of  the  vehicle.  These 
codes  are  very  costly  and  time  consuming  to  use.  While  these  approaches 
are  eminently  suited  to  the  final  "fine  tuning"  stage  of  the  vehicle 
design  process,  they  are  ill  suited  for  use  in  the  preliminary  design 
phase  where  a number  of  alternative  configurations  are  explored  and 
structural  details  remain  undefined.  The  present  elfort,  therefore, 
is  not  directed  toward  replacing  code  calculations  but  toward  supplementing 
them.  We  seek  a design-oriented  method,  rather  than  a capability  to 
analyze  fine  detail. 

At  present,  our  analysis  capability  is  restricted  to  axisymmetr ic 


response  of  projectiles  of  revolution  due  to  normal  impact  events.  Existing 
numerical  approaches  are  similarly  limited  because  three-dimensional 
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analysis  of  uncoupled  loads  is  not  yet  possible. 

In  addition  to  developing  the  new  approach,  it  has  been  applied 
to  the  penetrator  experiments  as  mentioned  earlier.  Reasonable  agree- 
ment between  predictions  and  measured  data  at  locations  on  the  primary 
penetrator  structure  is  an  important  benchmark  in  validating  the  new 
method  of  analysis. 

Basic  Approach 

The  impact  and  subsequent  earth  penetration  of  an  EPV  is  a complex 
type  of  phenomenon.  A simple  description  of  such  a phenomenon  is  only 
possible  if  drastic  assumptions  can  be  justifiably  made.  Several  are 
implicit  in  our  approach;  they  are  essent  al  to  achieving  our  objectives 
and  serve  as  the  groundwork  for  the  new  method  of  analysis. 

Considerable  effort  has  been  devoted  to  analyzing  the  impact  and 
penetration  of  earth  targets  under  the  assumption  that  the  vehicle  is 
rigid.  The  logic  here  is  that  the  loads  imposed  on  the  vehicle  by  the 
target  material  can  be  determined  in  this  way  as  a first  step.  This 
is  to  be  followed  by  an  analysis  of  the  vehicle  respose  to  these  known 
loads.  This  "uncoupling"  of  th.e  target  and  vehicle  response  obviously 
simplifies  the  design  analysis.  It  is  justifiable  if  the  vehicle  stiff- 
ness is  much  greater  than  that  of  the  target,  a condition  that  is  desirable, 
— if  not  necessary,  for  successful  penetration  to  great  depth.  Uncoupled 

load  determination  is  adopted  in  our  approach. 

A second  fundamental  assumption  that  underlies  the  practical  utility 
of  our  analysis  is  that  the  important  time  period  in  the  vehicle  response 
is  just  subsequent  to  impact.  We  choose,  therefore,  a stress  wave 
transmission  approach  to  analyze  vehicle  response.  Although  in  principle 
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a stress  wave  type  of  analysis  does  not  restrict  our  ability  to  predict 
long-time  response,  it  does  in  practice  as  following  reflections  from 
the  boundaries  of  the  vehicle  becomes  quite  cumbersome.  Furthermore, 
reflection  conditions  for  portions  of  the  vehicle  buried  in  the  target 
that  are  consistent  with  an  uncoupled  load  determination  approach  are 
subject  to  question  on  physical  grounds. 

Another  important  assumption  is  that  the  essential  nature  of  the 
vehicle  response  is  longitudinal  and  simple  bar-like  due  to  normal  im- 
pact. This  implies  that  one-dimensional  plane  waves  only  need  be  con 
sidered  and  that  the  stress  state  is  uniaxial.  Rotational  symmetry 
of  the  vehicle  about  a longitudinal  axis  and  a degree  of  slenderness 
of  the  penetrator  are  necessary  for  this  assumption  to  be  valid,  in 
addition  to  a perfectly  normal  impact/penetration  event.  A further 
simplification  results  if  only  Hookean  elastic  material  response  is 
considered;  this  is  consistent  with  the  use  of  very  brittle,  high  strength 
steels  for  high  velocity  penetrator  bodies. 

The  analysis  is  conducted  within  the  general  framework  defined 

by  the  above  assumptions  and  restrictions.  The  approach,  although 

general,  requires  a concrete  example,  such  as  that  provided  by  the 

3 

penetrator  experiments  conducted  by  Avco  , to  convey  the  role  of  judg 
ment  in  obtaining  good  results. 

METHOD  OF  ANALYSIS 

Preliminary  Remarks 

There  are  three  modeling  tasks  to  be  completed  in  order  to  charac- 
terize the  impac t/penetrat ion  event  and  the  EPV.  The  first  is  defining 
the  loading  environment  to  which  the  EPV  is  exposed;  this  must  be  done 
in  a manner  that  is  consistent  with  the  simplicity  desired  and  the  one- 
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dimensional  bar-like  representation  of  the  penetrator  primary  structure. 
An  obvious  choice  is  to  consider  the  resultant  axial  force  exerted  by 


I 

■ 
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the  target  as  acting  at  a representative  station  on  the  axis  of  the  EPV. 

Secondly,  as  indicated  by  the  sketch  in  Figure  1,  most  EPV's  can  be 
viewed  as  being  composed  of  two  distinct  portions.  There  is  a forebody 
which  is  highly  tapered  geometrically*that  includes  the  forward  20  per 
cent  (approximately)  of  the  total  length.  The  remaining  portion  of 
the  vehicle,  the  aftbody,  is  rather  uniform  by  comparison.  It  is  natural, 
therefore,  to  model  these  two  regions  in  different  ways  appropriate 
to  their  respective  physical  characteristics.  Consequently,  idealized 
forebody  and  aftbody  models  are  defined  and  analyzed  separately,  as 
will  be  shown  in  the  sequel. 

Judgment  on  the  part  of  the  analyst  is  exercised  in  defining  the 

loading,  forebody  and  aftbody  models.  This  is  always  the  situation 

if  good  results  are  to  be  obtained  by  a simple  method  of  analysis. 

In  the  sequel,  we  present  results  that  are  the  consequences  of  our 

3 

judgement  in  describing  the  penetrator  experiments.  It  is  beyond 
the  scope  of  our  work  to  consider  the  sensitivity  of  predictions  to 
modeling  details. 

Loading  Idealization 

The  uncoupled  loads  on  the  penetrator  for  impact  into  Dakota  sand- 

4 

stone  at  1800  feet  per  second  were  calculated.  A two-dimensional 
finite  difference  Lagrangian  computer  code  was  used.  Results  are  pre- 
sented in  Reference  4 for  pressure  distributions  at  various  times,  time 
histories  of  pressure  at  fixed  locations  and  several  overall  parameters 
*Some  EPV's  are  tapered  over  their  entire  length  (all  forebody). 
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as  functions  of  time.  Figure  6 of  Reference  4 is  of  paramount 

importance  to  the  present  analysis;  it  shows  the  calculated  resultant 
axial  force  time  history  which  serves  as  the  input  to  our  two  structural 
mode  Is . 

The  axial  force  time  history  in  Figure  6 of  Reference  4 contains 
the  typical  "noise"  or  oscillations  characteristic  of  the  calculation 
method.  For  our  purposes,  the  time  history  is  represented  by  a series 
of  line  segments  as  shown  schematically  in  Figure  2.  The  force  drops 
after  time  t in  Figure  2 because  the  target  blocks  used  in  the  tests, 
being  of  finite  extent,  have  a free  outer  surface;  a relief  wave,  there- 
fore, is  created  by  a free  boundary  reflection. 

A study  of  Figure  6 of  Reference  4 indicates  that  our  idealized 
line  segment  time  history  model  will  correspond  well  if  we  choose 

F = 43,000  lbs.  t,  = 160  microseconds 

o 1 

F^  = 314,000  lbs.  t^  = 205  microseconds 

F^  = 47,300  lbs.  t^  = 400  microseconds 

(The  above  force  and  time  parameters  are  defined  in  Figure  2.)  It  takes 
a finite  time  after  impact  for  the  axial  force  to  develop  appreciable 
- magnitude.  Figure  6 of  Reference  4 indicates  that  a 10-microsecond 

delay  in  response  is  appropriate  to  account  for  this  effect  as  F does 
not  correspond  to  zero  time  (impact). 

It  is  also  necessary  to  establish  the  point  of  application  for 
the  resultant  force  on  the  penetrator's  longitudinal  axis  of  symmetry. 

In  reality,  this  "center  of  force"  is  time  dependent.  It  moves  from 
the  tip  of  the  penetrator  nose  at  the  instant  of  impact  to  its1  aftmost 
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RESULTANT  AXIAL  FORCE 


t,  TIME 

Figure  2.  Idealized  Force  Time  History 


10 


i 


* 


i 


' 

I 


location  (approximately  station  2.05)  about  250  microseconds  later^ 
as  the  penetrator  is  progressively  encapsulated  by  the  target  material. 
Since  accounting  for  this  shift  of  the  center  of  force  with  time  com- 
plicates the  analysis  beyond  the  intended  scope,  a simplification  is 
introduced.  The  force  is  assumed  to  always  act  at  a fixed  location. 

This  assumption  implies  that  response  variables  corresponding  to  points 
forward  of  the  aftmost  position  of  the  center  of  force  cannot  be  estimated 
on  a rational  basis.  In  view  of  this,  the  aftmost  location,  station 
2.05,  is  selected  as  the  fixed  location. 

The  above  approach,  although  there  is  no  uniqueness  associated 
with  our  choice  for  the  effective  center  of  force,  is  the  simplest  way 
to  treat  the  progressive  encapsulation  by  the  target  material.  Any 
other  approach  involves  time-wise  numerical  solution  of  equations,  which 
is  certainly  to  be  avoided  if  the  original  objectives  are  to  be  met. 
Detailed  information  on  the  nose  tip  region  is  sacrificed  in  order  to 
obtain  the  great  simplicity  demonstrated  in  the  following  sections. 


Forebody  Model  and  Analysis 

The  forebody  portion  of  the  EPV  is  characterized  by  a high  degree 
of  geometric  taper.  Ogive  shapes  are  the  usual  choice.  A sketch  of 
our  forebody  model  is  shown  in  Figure  3.  The  origin  for  the  longitu- 
dinal coordinate  x (positive  aftward)  is  taken  at  the  assumed  center 
of  force.  Note  that  the  nose  portion  forward  of  this  point  is  not  load- 
ed according  to  this  idealized  model;  it  contributes,  therefore,  only 
a mass  M which  acts  effectively  at  the  origin.  Uniaxial,  simple  bar- 
like response  is  characterized  by  the  single  longitudinal  displacement 
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u(x,t)  which  describes  the  translation  of  the  cross  section  at  x for 
any  time  t. 

The  differential  equation  governing  longitudinal  motion  of  a homo- 
geneous bar  with  varying  cross  sectional  area  A is 

( EAu,  ),  - pAu,  =0  (1) 

x x ’tt 

The  elastic  modulus  is  E and  p is  the  mass  density  of  the  bar  material. 
This  equation  is  valid  provided  that  the  bar  is  slender  and  the  variation 
of  A is  smooth  and  not  too  severe;  if  these  conditions  are  not  met, 
three-dimensional  effects  become  important,  thereby  tending  to  undermine 
the  accuracy  of  this  simple  theory. 

A further  simplification  results  if  we  restrict  out  attention  to 
tapered  bars  of  revolution  with  a linearly  varying  outer  radius;  in 
this  case,  an  exact  solution  has  been  found  to  Equation  (1).  This 
"conical"  approximation  describes  the  major,  first  order  influence  of 
geometric  taper.  The  radius  r(x),  then,  satisfies  an  equation  of  the 
form 


r(x)  = a + a x (2) 

The  radius  at  the  origin  is  a,  and  a is  a parameter  characterizing 
the  taper.  The  cross  sectional  area  is 


A(x)  = n (a  + a x)‘ 


Equation  (1)  can  be  written  in  the  following  form: 


xx 


2«  1 _ n 

(a  + ax)  u,x  2 u,tt 
c 


(3) 


(4) 
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The  velocity  of  propagation  of  di 
The  solution  to  Equation  (4) 
x-direction  is 


sturbances  in  the  bar  is  c = \/e7"p. 
for  a wave  traveling  in  the  positive 


u(x, t ) 


a f(t  - -) 

(a  + » x) 


(5) 


f is  an  arbitrary  function  to  be  determined  by  the  boundary  condition 
at  x = 0.  The  validity  of  this  solution  is  easily  verified  by  direct 
substitution  into  Equation  (4). 

The  appropriate  boundary  condition  to  be  applied  to  x = 0 corresponds 
to  the  equation  of  motion  for  the  mass  M.  It  may  be  written  as  follows 

M u, tt (0, t ) - F(t)  ♦ EA(0)  u,x(0,t)  (6> 

With  the  aid  of  Equations  (3)  and  (3),  this  result  may  be  simplified 
to 

2 

Mf"  (t)  ♦ f'(t)  ♦ tt  a aEf(t)  = F(t)  (7) 

c 


The  notation  ( )'  is  used  to  indicate  an  ordinary  derivative  with  respect 
to  the  argument  of  the  function.  The  mass  M is  small  compared  to  the 
total  mass  of  the  model.  Consequently,  a simplification  results  if 
the  M-term  is  neglected  in  (7).  The  following  simplified  equation 
results : 


f'(t)  + <JL£  f(t) 
a 


F(t) 


(8) 
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This  equation  reflects  the  effect  of  taper  by  the  presence  of  the  para- 
meter a 


The  solution  to  Equation  (8)  can  be  written  for  arbitrary  F(t) 
with  the  aid  of  the  indicial  admittance  function  A(t).  A(t)  is  defined 
as  the  solution  corresponding  to  a unit  step  function  force  distribution; 
it  is  * 


A (t) 


1 

HaaE 


(1 


(9) 


The  general  solution  for  arbitrary  F(t)  is 

t 

f(t)  = J F(T)  A'(t-T)  dT  (10) 

0 


For  early  times  (ts  t^),  the  force  distribution  according  to  our 
loading  model  shown  in  Figure  2 can  be  written  as 


F(t) 


(F, 


F + 
0 


- V 


(t  s tL) 


For  this  distribution 


f(t) 


ttaaE 


t 


) 


r 


+ 


(F1  - Fo> 


tjJTor  cE 


«c 

1 + e a 


t 


(ID 


(12) 


(12) 


The  axial  strain  « (x,t)  can  be  determined  from  f(t  - — ) and  Equation 
x c M 

(5).  The  result  can  be  written  as  follows: 

★ ...  w . 

The  initial  condition  A(0)  = 0 has  been  used  to  obtain  this  result. 
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(13) 


(x , t ) = u , 

X X 


= Cje 


(t  - i) 


+ C„ 


c,(t  - -) 

3 c 


The  functions  Cj  , and  are 


II 

o 

(i  - 

1 ) 

l TTE  ( a + a x ) 

\ a 

(a  + ax) I 

1 

1 

© 

a 

1 

t c j . 

1 (a  + 

. 2 
ax) 

(a  + ax) 

0 

f E(a  + «x)" 


<fi-v  r i 

t*E3  c Ta  ~ 


a x)  , _ ,2 

(a  + Qfx) 


(14) 


(13) 


C3 


<F1  - Fo> 

t nE(a  + ax)' 


(16) 


In  applying  the  above  results  to  the  penetrator  shown  in  Figure 
1,  the  fact  that  the  nose  region  is  composed  of  two  materials  has  been 
ignored.  The  properties  of  the  steel  outer  shell  have  been  used.  The 
appropriate  parameter  values  are 

a = 0.2126 

c = 197,174  inches/sec 
a = 0.874  inches 
E = 28.5  x 10^  psi 

It  will  be  apparent  from  the  discussion  of  results  that  this  approxi- 
mation for  the  forebody  provides  excellent  strain  predictions. 
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Aftbody  Model  and  Analysis 

In  contrast  to  the  EPV  forebody,  the  primary  load  carrying  structure 
of  the  aftbody  is  relatively  uniform  in  geometry  and  stiffness.  As 
a first  approximation,  therefore,  the  analytical  model  is  based  upon 
uniform  or  nominal  stiffness  for  the  entire  length  of  the  aftbody. 

The  forebody's  influence  is  confined  to  introducing  a (usually  substan- 
tial) mass  M in  the  most  forward  location  and  a path  for  transmitting 
the  resultant  force  F.  This  idealized  model  is  shown  in  Figure  4. 

There  are  two  uniform  elastic  bars  which  act  in  parallel  connected  to 
M . This  is  because  the  outer  body  shell  and  inner  liner  or  payload 
shell  of  the  penetrator  sh  wn  in  Figure  1 are  separated  by  an  air  gap 
and  do  act  as  independent  load  transmission  paths. 

We  use  the  subscript  "1"  for  the  bar  corresponding  to  the  outer 

body  shell  and  "2"  for  the  payload  shell.  The  structural  effect  of 

the  potting  material  used  in  the  payload  module  (Figure  1)  is  neglected 

completely.  Since  the  bars  are  assumed  to  be  uniform,  traveling  waves 

in  the  positive  x-direction  in  the  bars  are  represented  by  the  functions 

g (t  - — ) and  g (t  - — ).  Physically  both  bars  are  attached  to  the 
1 C1  2 c 2 

mass  at  the  origin;  therefore 

gL ( t ) = g2(t)  = g(t)  (17) 

A boundary  condition,  similar  to  Equation  (6)  and  expressing  the 
equation  of  motion  for  M , can  be  formulated  for  use  in  determining 
g.  It  is 


M g"(t)  + 

o 


g'(t)  = F(t) 


1 7 


(18) 


Figure  4.  Aftbody  Model 
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Also,  we  require 


g(0)  = g'(0)  = 0 


(19) 


as  the  bars  are  at  rest  initially. 

It  is  convenient  to  introduce  the  parameter  p defined  as 


P = 


(20) 


Equation  (18),  then,  is  equivalent  to 


g"  + Pg'  = 


F(t) 


(21) 


By  analogy  with  the  solution  process  used  to  determine  f for  the  fore- 
body model,  we  seek  the  response  due  to  a unit  step  function  force. 

Let  v be  the  velocity  (g1)  for  a unit  step  force  input;  it  is  given 
by  the  following  expression: 


v 


1 


( 1 - e 


pt . 


(22) 


The  velocity  corresponding  to  an  arbitrary  force  F(t)  is,  therefore 

v =g'  = f F(T  ) v'(t  -T)dT  = i-  f e_M,(t  ~T)F(t)  dT  (23) 

J M J 

0 0 

The  force  time  history  model  defined  in  Figure  2 is  characterized 
by  three  distinct  analytical  expressions,  one  corresponding  to  each 
line  segment.  They  are 
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1 


F(t)  = 

F + 
0 

(F1  - f 

oL, 

(0  ^ t * t ) 

(ID 

ei 

F ( t ) 

" F1 

s t s t2) 

(24) 

F(t)  = 

(F1 

" V 

(t  - 

t2)  (t2‘tst3) 

(25) 

Tt^ 

- ‘2> 

These  expressions  may 

be  substituted 

into 

(23)  to  determine  the 

particle 

velocity  v.  The  axial 

strains,  then 

, in 

bars  1 and  2 are  found 

from 

the  following  equations: 

€ 

xl 

8l,x  = 

_ 1 

C2 

v(t  - - ) 

C1 

(26) 

e 

x2 

g2,x  = 

1 

’ "2 

v(t  - - ) 
C2 

(27) 

These  theoretical  results  have  been  applied  to  the  aftbody  (the 
portion  of  the  penetrator  aft  of  station  8)  of  the  penetrator  shown 
in  Figure  1.  The  following  parameter  values  have  been  used  in  the 
ana  lysis: 

= 26.5  x 10^  psi 
c 2 = 188 , 141  in/ sec 
= 0.4  in^ 
p.  = 11,871  sec"1 

The  results  for  a semi -infinite  aftbody  are* 

*The  units  of  velocity  used  are  inches  per  second. 


Ej  = 28.5  x 10  psi 


c | = 197,174  in/ sec 


A = 3.0  in 
Mq  = 0.0415 


lb . -sec 
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V - 87.4  ( 1 - e C)  + 289.7  ut 

L 

- (1  - 

J 

(0  1 1 ) 

(28) 

v = e""*  L 2533.9  + 638.1  (e*1 

- 6.7) 

3 (t. 

£ t s t2) 

(29) 

v = 5541.0  e~“l  +638.1  (1  - e”k( 1 ~ t2 * 

- 25.9 

, ( - 

1 ) - 1.4  e*  u 

(t  - 

(t2 

‘ C3) 

(30) 

A graphical  representation  of  the  strain  in  bar  1 obtained  using  Equations 
(28)  - (30)  appears  in  Figure  5. 

As  is  usual  in  wave  transmission  analyses  of  the  present  kind, 
baseline  information  is  provided  by  the  semi -in f in i * e mode!  solution 
described  above  and  shown  in  Figure  5.  The  effect  of  finite  penetrator 
length  is  accounted  for  by  appropriate  reflection  conditions  at  the 
aft  end.  Superposition  of  traveling  wave  solutions  is  used  to  satisfy 
the  aft  end  boundary  condition,  thereby  providing  an  exact  solution 
for  dynamic  response  consistent  with  the  loading  and  penetrator  models. 

A convenient  method  of  treating  reflections  is  described  in  the  Appendix. 

Synopsis  of  the  Method 

Our  method  of  analysis  is  summarized  briefly  as  follows: 

(1)  Define  loading,  forebody  and  aftbody  models; 

(2)  Determine  the  response  of  the  forebody  and  aftbody  models 
as  if  they  are  semi -in f ini te  in  length; 

(3)  Use  superposition  of  the  semi - i n f in i te  traveling  wave  solu- 
tions to  satisfy  appropriate  boundary  conditions  that  correspond 
to  finite  penetrator  Length. 

Corresponding  to  any  model  definitions,  the  response  is  obtained  from 
simple  closed  form  expressions  and  their  superposition.  As  will  be 
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shown  next,  this  simplicity  is  also  accompanied  by  excellent  agreement 
with  available  experimental  data  on  response  of  the  primary  structure 


RESULTS  AND  DISCUSSION 

Preliminary  Remarks 

From  the  outset,  our  primary  concern  has  been  predicting  the  response 
of  the  load  carrying  structural  elements  by  simple  means.  Simplicity 
and  ease  of  application  have  certainly  been  achieved.  We  now  show  that 
our  approach  also  provides  excellent  quantitative  predictions  for  strains 
and  accelerations. 

In  the  following,  strain  and  acceleration  predictions  are  discussed 
separately,  as  are  the  response  of  the  penetrator  main  body  outer  shell 
and  the  internal  structure.  This  way  of  categorizing  the  results  is 
both  convenient  and  useful  in  judging  the  validity  of  the  modeling 
assumptions . 

The  analytical  results  presented  are  directly  comparable  to  the 
experimental  ones  that  are  given  in  Reference  3.  While  our  intent  is 
not  to  provide  an  exhaustive  correlation  study,  representative  experi- 
mental data  are  presented  on  the  same  figures  with  theoretical  predictions 
in  several  instances.  These  selected  benchmarks  are  sufficient  to 
illustrate  the  generally  fine  degree  of  agreement  between  theory  and 
experiment  that  has  been  achieved. 

Outer  Shell  Strain  Predictions 

Axial  strains  have  been  measured  at  four  locations  on  the  exterior 
of  the  penetrator  main  body  outer  shells  during  the  two  tests. 


Station  4, 


the  most  forward  location,  is  the  only  site  on  the  penetrator  forebody, 
lhe  remaining  sites  are  stations  9,  19.5  and  2 7.3  on  the  aft  body, 
lheoretical  axial  strain  time  histories  determined  by  our  analysis 
method  for  these  locations  appear  in  Figures  6-9. 

The  strain  gages  at  station  4 are  wiped  off  during  tie  penetrat 
of  the  target.  Consequently,  only  short-time  data  is  available, 
this  case,  only  the  force  representation  in  Equation  (11)  n«  . d •>. 
sidered.  From  a study  of  Figure  1,  it  is  to  be  expected  that  r<  • • t 
of  the  primary  stress  wave  will  be  produced  at  station  8 where  the  .r  , 
ballast  plug  ends.  Station  8 behaves  approximately  as  a tree  end  1 r 
the  forebody,  while  the  buried  portion  of  the  penetrator  nose  can  be 
considered  fixed.  As  a consequence  of  these  double  reflections,  the 
experimental  data  in  Figure  6 exhibits  two  noticeable  pauses  in  the 
otherwise  monotonic  trend.  The  theoretical  prediction  is  for  a semi- 
infinite forebody,  so  the  pauses  do  not  occur  in  that  curve.  Clearly, 
there  is  excellent  agreement  otherwise,  and  it  is  apparent  that  the 
first  passage  semi-infinite  prediction  is  a conservative  design  esti- 
mate . 

Response  at  station  4 has  been  studied  extensively.  The  theoret- 
ical results  presented  in  Figure  6 account  for  the  mass  M (Equation 
(7)).  Another  solution  based  upon  Equations  (8)  and  (13)  is  in  close 
agreement.  Attempts  to  account  for  the  double  reflection  pauses  have 
met  with  reasonable  success  also.  For  brevity,  only  the  one  prediction 
i s shown . 

The  strain  predictions  at  stations  9,  19.5  and  27.3  in  Figures 
7-9  have  all  been  obtained  from  the  semi-infinite  aftbody  model  response 
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shown  in  Figure  5.  The  aft  end  of  the  outer  shell  at  station  32  behaves 
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essentially  as  a free  end.  Consequently,  a relief  wave  is  generated 

which  is  opposite  in  sign  (tensile)  to  the  incident  (compressive)  wave; 

a simple  graphical  approach  to  superimposing  the  effects  of  these  two 

waves  is  given  in  the  Appendix.  The  arrival  >f  the  relief  wave  at  each 

station  is  easily  identified  in  the  figures. 

The  prediction  for  station  19.5  is  compared  with  experimental  data 

in  Figure  8.  Note  that  a shift  to  the  left  of  25  microseconds  in  the 

theoretical  curve  would  produce  nearly  perfect  agreement  between  theory 

and  experiment.  The  time  scale  of  the  experimental  data  is  only  accurate 

3 

to  within  this  time  increment.  Although  no  effort  has  been  made  to 
adjust  the  time  origin  for  improved  agreement,  this  particular  case 
is  conspicuous  in  that  a shift  produces  striking  agreement. 

Strain  Predictions  For  Internal  Structure 

The  payload  shell  or  liner  is  the  primary  load  carrying  element 
in  the  interior  of  the  penetrator.  The  shell  is  not  continuous,  but 
segmented  in  a modular  fashion.  The  segments  are  threaded  at  the  ends 
with  the  joints  reinforced  by  rings. * It  is  to  be  expected,  therefore, 
that  variability  in  these  mechanical  joints  makes  the  precise  nature 
of  the  load  transmission  path  indeterminant  in  practice. 

Our  aftbody  model  is  based  upon  continuity  of  the  payload  shell; 
the  joints  and  small  reinforcing  rings  are  totally  ignored.  There  are 
two  very  large  titanium  rings*  located  at  stations  23.25  and  31.25 
(Figure  1)  which  cannot  be  ignored,  however. 

The  two  titanium  rings  have  been  modeled  as  concentrated  masses 
attached  to  the  payload  shell.  The  addition  of  a concentrated  mass 
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to  an  elastic  bar  is  analyzed  in  detail  in  Reference  5.  Initially, 
the  concentrated  mass  produces  a reflection  of  the  incident  stress  wave 
as  if  the  mass  is  a fixed  end;  the  reflection  decays  exponentially  with 
time  thereafter.  Simultaneously,  the  stress  wave  transmitted  beyond 
the  mass  builds  up  exponentially  until  the  original  incident  stress 
level  is  reached.  For  stations  forward  of  the  first  ring,  it  is  con- 
servative to  assume  that  a fixed  end  reflection  occurs. 

Axial  strain  data  has  been  obtained  at  two  stations  on  the  inside 
of  the  payload  shell,  stations  9 and  19.5,  forward  of  the  first  titanium 
ring.  The  theoretical  predictions  shown  in  Figures  10  and  11  are  based 
upon  the  assumption  of  a fixed  end  at  the  first  ring.  The  station  9 
prediction,  although  slightly  conservative,  agrees  very  well  with  the 
experimental  data  from  Reference  3.  In  Figure  11,  representative  ex 
perimental  data  is  also  presented  for  station  19.5;  note  that  the  pre- 
diction indicates  a more  rapid  strain  build-up  than  the  measured  data 
(a  conservative  trend),  although  the  peak  strain  values  are  in  good 
agreement . 

Axial  strain  has  been  measured  also  at  station  31.5  on  the  inside 

of  the  payload  shell.  This  location  is  aft  of  both  large  titanium 

rings.  The  measured  strain  is  much  smaller  in  magnitude  than  the  others 

3 

and  opposite  in  sign  (tensile)  . Our  modeling  approach  to  these  rings 
would  suggest  that  the  strain  at  this  location  is  essentially  zero  in 
the  time  period  of  interest;  consequently,  there  is  no  predicted  time 
history  included. 

In  addition  to  the  payload  shell,  strain  gages  were  mounted  on 
the  flexible  potting  material  at  Station  16.  Our  aftbody  model  does 
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not  account  for  any  contribution  by  the  potting  to  structural  stiffness 
as  it  is  quite  small.  The  axial  strain  in  the  potting  can  be  estimated, 
however,  if  strain  compatibility  with  the  payload  shell  is  assumed. 

This  approximation  results  in  a potting  axial  strain  at  Station  16  that 
is  essentially  identical  to  that  in  the  payload  shell  at  station  19.5 
with  a slight  shift  of  the  time  axis;  the  same  prediction  curve  suffices, 
therefore.  This  is  in  good  agreement  with  the  experimental  results 
given  in  Reference  3. 

Hoop  Strain  Predictions 

All  hoop  strains  corre'pond  to  simple  Poisson  expansion  only  accord- 
ing to  our  simple  bar-like  nodels.  Consequently,  separate  results  for 

the  hoop  strains  are  not  presented.  This  simple  uniaxial  stress  state 

3 

is  a good  approximation  over  most  of  the  penetrator.  It  cannot  be 
applied  at  Station  4,  however,  where  three-dimensional  effects  due  to 
nose  encapsulation  obviously  occur;  prediction  of  these  effects  is 
beyond  the  scope  of  our  approach. 

Acceleration  Predictions 

Two  accelerometers  were  mounted  on  primary  load  carrying  structure 
for  the  Avco  tests,  one  forward  at  Station  7 and  one  aft  at  station 
32  . Acceleration  predictions  for  these  locations  are  presented  in 
Figures  12  and  13,  respectively  . These  predictions  have  been  determined 
using  the  aftbody  model. 

Also  shown  in  Figure  13  is  data  from  Reference  3 for  comparison. 

The  trends  correspond  quite  well,  but  a high  frequency  oscillation  in 
the  test  data  is  not  predicted  by  our  analysis. 
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CONCLUSIONS  AND  RECOMMENDATIONS 

The  work  presented  here  is  a pioneering  effort  to  develop  a simple, 
reliable  preliminary  design  analysis  method  for  normal  impact  of  eartt 
penetrators.  In  addition,  the  method  has  been  applied  to  predict  the 
outcome  of  actual  impact/penetration  tests.  The  results  obtained  for 
this  particular  application  are  in  excellent  agreement  with  available 
experimental  data  associated  with  locations  on  the  primary  load  carrying 
structure  of  the  penetrators  tested.  We  conclude,  therefore,  that  the 
objectives  stated  in  the  introduction  have  been  met.  We  believe  that 
the  approach  is  basically  sound  and  that  the  evidence  supporting  this 
conclusion,  although  limited  to  only  one  penetrator  configuration,  is 
never-the-less  convincing. 

The  method,  although  useful  as  far  as  it  goes,  is  still  not  yet 
beyond  the  embryonic  stage.  Additional  work  needs  to  be  done,  and 
several  questions  remain  unanswered.  The  following  items  are  examples 
of  areas  that  require  attention: 

(1)  Explore  alternatives  to  the  fixed  center  of  force  approach 
which  better  describe  the  dynamic  encapsulation  of  the  pene- 
trator by  the  target  material; 

(2)  Study  the  limitations  on  the  applicability  of  the  one-dimen- 
sional models;  (e.g.,  how  slender  must  the  penetrator  be  for 
the  bar-like  models  to  be  valid?) 

(3)  Develop  simple  loading  prediction  methods  which  do  not  require 
large-scale  numerical  calculations  or  which  are  not  empirical 
in  nature.  These  methods  should  be  in  the  same  spirit  as 

the  present  penetrator  response  analysis; 

(4)  Verify  the  validity  of  the  method  for  penetrators  which  differ 
substantially  in  configuration  from  the  ones  considered  here; 
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(5)  Refine  the  forebody  model  to  account  for  second  order  effects 
of  taper; 

(6)  Develop  sim'lar  models  and  an  analysis  method  that  applies 
to  non-normal  impact  events; 

(7)  Develop  solutions  valid  for  long-time  response  and  match  them 
to  the  short-term  wave  analysis. 

The  above  listing  order  does  not  imply  an  order  of  priority.  Successful 
accomplishment  of  these  tasks  will  insure  that  the  desired  design  analysis 
capability  will  be  effective  and  accepted  by  penetrator  analysts  and 
designers . 
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APPENDIX 


A Simple  Graphical  Approach  for  Wave  Reflection  Analysis 

As  indicated  in  the  discussion  of  the  aftbody  model,  baseline 
information  for  a wave  transmission  analysis  is  provided  by  a semi- 
infinite model  solution  as  shown  in  Figure  5.  The  effect  of  finite 
penetrator  length  is  accounted  for  by  appropriate  reflection  conditions. 
For  clarity,  we  will  outline  the  approach  used  to  satisfy  free  end 
conditions  at  the  aft  end  of  the  main  body  outer  shell. 

At  a free  end,  the  axial  stress,  and  hence  the  axial  strain,  is 
zero  for  all  time.  This  condition  is  satisfied  if  a wave  of  opposite 
sign  but  same  geometric  form  as  the  incident  wave  travels  in  the  opposite 
direction  in  the  bar;  the  incident  semi-infinite  model  wave  and  the 
"reflected"  wave  of  opposite  sign  must  be  timed  such  that  the  sum  of 
the  axial  strains  due  to  both  waves  is  always  zero  at  the  end  point. 

The  strain  distribution  within  the  penetrator  at  any  time  prior  to  the 
time  that  the  reflected  wave  reaches  the  forward  end  is  the  sum  of  the 
contributions  due  to  the  two  waves. 

This  superposition  is  easily  accomplished  graphically  as  follows. 

For  a given  location  on  the  penetrator,  the  arrival  times  for  the  in 
cident  semi -in f in i te  solution  strain  wave  and  the  free  end  reflected 
wave  of  opposite  sign  are  determined  for  the  wave  speed  of  the  material. 
Between  these  two  times,  the  solution  is  the  same  form  as  the  semi- 
infinite solution  with  origin  at  the  former  time.  For  times  after  the 
latter  arrival  time,  appropriate  increments  must  be  subtracted  from 
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the  incident  wave  strain  values  to  correspond  to  a reflected  wave  moving 
in  the  opposite  direction;  the  origin  of  the  reflected  wave  is  the 
latter  arrival  time.  This  superposition  is  easily  accomplished  with 
tlie  aid  of  a light  table,  copies  of  the  semi -i n f in i te  incident  wave 
profile,  and  dividers.  The  resulting  profile  for  a particular  station 
can  be  determined  quite  rapidly  in  this  way. 

Fixed  end  conditions  require  addition  of  like  strain  wave  profiles, 
rather  than  subtraction.  Acceleration  profiles  may  be  obtained  in  a 
similar  manner,  but  the  role  of  addition  and  subtraction  are  exactly 
reversed  for  accelerations  due  to  free  and  fixed  ends. 
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